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We propose that repeated photoexcitation/ionization of high Z atoms of highly 

relativistic flows by star-light in dense stellar regions followed by emission of de- 
' cay /recombination photons, which are beamed and boosted to 7 ray energies in 
the observer frame, produce gamma ray bursts (GRBs). We show that this over- 
looked mechanism is able to convert efficiently baryonic kinetic energy release 
in merger or accretion induced collapse of neutron stars into cosmological GRBs 
and reproduces remarkably well all the main observed properties of GRBs. 

The origin of gamma ray bursts (GRBs), which have been discovered 35 years ago, is still 
a complete mystery [1]. Their observed isotropy in the sky, deficiency of faint bursts and the 
lack of concentration towards the Galactic center, in the Galactic disk and in the direction 
of M31, strongly suggest [2] that they are cosmological in origin [3]. A cosmological origin 
implies [4] that GRBs have enormous luminosities during short periods of time, 
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where d = cfoslO 28 cm is their luminosity distance, <p — 10~ 6 </>6 ergs cm~ 2 s~ 1 is their measured 
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energy flux, x 1S a bandwidth correction factor and Af2 is the solid angle into which their 
emission is beamed. Moreover, the short durations of GRBs imply very compact sources. 



2 ' Relativistic beaming is required then, both in order to explain the absence of accompanying 
optical light and X-ray emission, and in order to avoid self opaqueness due to 77 — > e + e~ . All 
these considerations [4] seem to support the favorite cosmological model of GRBs; relativistic 
fireballs [5] formed by mergers of neutron stars (NS) or neutron stars and black holes (BH) 
in close binary systems [6] due to gravitational wave emission [7], or by accretion induced 
collapse (AIC) of NS and white dwarfs (WD) [6,8]. Indeed, the above considerations and 
the observed rate of GRBs favor NS mergers/AIC as the source of cosmological GRBs. 
Nevertheless, no mechanism has been convincingly shown to be able to convert a large 
enough fraction of their binding energy release into 7 rays and to explain simultaneously the 
observed complex light curves, duration distribution and spectral behavior of GRBs [4]. 

Most of the binding energy released in NS mergers/AIC is expected to be in the form of 
neutrinos, gravitational waves and kinetic energy of ejected material [9]. No mechanism is 
known which converts efficiently gravitational waves into 7 rays. Neutrino annihilation [6] 
[yv — > e + e~), and neutrino pair production in strong magnetic fields (wy v — > ve + e~) near 
merging/collapsing NS cannot convert enough binding energy into e + e~7 fireballs or rela- 
tivistic e ± beams which can produce cosmological GRBs, because of baryon contamination 
[10] and the low efficiency of these processes [11]. Thus, if NS-NS and NS-BH mergers, or 
AIC of WD and NS, produce GRBs, the production must proceed through conversion of 
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baryonic kinetic energy of a highly relativistic flow (a fireball or a jet) into 7 rays. In this 
letter we propose that repeated photoexcitation/ionization of the highly relativistic atoms 
of the flow by star-light in dense stellar regions followed by emission of decay/recombination 
photons which are beamed and boosted to 7 ray energies in the observer frame (see Fig. 1), 
produce cosmological GRBs. We show that this simple mechanism, which has been over- 
looked, is able to convert enough baryonic kinetic energy of highly relativistic flows in dense 
stellar regions into cosmological GRBs. We also show that it predicts remarkably well the 
main observed properties of GRBs; their burst size, duration distribution, complex light 
curves and spectral evolution. 

The surface of neutron stars is believed to consist of iron-like nuclei [9]. Therefore, 
we assume that the ejected mass in NS mergers or in AIC of NS or WD contains such 
high Z nuclei. In view of the uncertainties in modeling merger/AIC of compact stellar 
objects [11], rather than relying on numerical simulations, we deduce the total relativistic 
kinetic energy release in such events from observations of type II supernova explosions, which 
are driven by gravitational core collapse into NS or BH. In type II supernova explosions, 
typically, 1OM are accelerated to a final velocity of v ~ 10 4 km s' 1 [9], i.e., to a total 
final momentum P ~ 10M o t>. Since core collapse is not affected directly by the surrounding 
stellar envelope, we assume that in merger/AIC collapse of NS or WD a similar impulse, 
/ Fdt PS P, is imparted to the ejected mass. If the ejected mass is much smaller than a 
solar mass, AM <C M , then it is accelerated to a highly relativistic velocity and its kinetic 
energy is given approximately by E K ~ Pc ~ 6 x 10 53 erg. More than 1CT 3 of this highly 
relativistic baryonic kinetic energy must be converted into 7 rays in order to produce a 
typical cosmological GRB of ~ 10 51 erg. 

The natural birth places of close binary systems are the very dense stellar regions in 
galactic cores and collapsed cores of globular clusters (GC). These GC regions have very large 
photon densities and column densities, n 7 ~ 3L/47icR 2 e and iV 7 ~ L/ncRe , respectively. 
For instance, the core of our Milky Way (MW) galaxy has a surface brightness, S ~ 1.2 x 
10 7 L Q pc~ 2 , typical photon energies, e ~ 1 eV, and a radius, R ~ 0.5 pc [12]. These values 
yield, on average, iV 7 ~ 10 23 cm' 2 . (Actually, this value underestimates the column density 
because we have neglected the contribution from stars beyond R. Its inclusion yields iV 7 that 
are several times larger. A similar value is obtained for the nearby M31 galaxy. Much larger 
values are obtained for active galactic nuclei). The highly relativistic ejecta expands, cools 
and becomes transparent to its own radiation after a relatively short time [5]. Its partially 
ionized, highly relativistic atoms see the interstellar photons blue shifted to X-ray energies. 
These photons photoexcite/photoionize the high Z atoms of the flow which subsequently 
decay /recombine radiatively by isotropic X-ray emission in the flow's rest frame. For a 
highly relativistic flow with a Lorentz factor T ^> 1, the isotropic X-rays are beamed along 
the flow, and as long as 6 < 1/T their energies, ex, are boosted to 7 ray energies, 

e^2Te x /(l + T 2 6 2 ), (2) 
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where 9 is the angle of the flow direction relative to the observer. The differential fraction 
of the emission that is directed towards the observer is given by 



dI/d(cos6) « 2T 2 /(1 + 2T 2 6 2 ). (3) 

The total baryonic kinetic energy which is converted into beamed 7 ray emission by 
repeated excitation/decay can be estimated as follows: The total photoabsorption cross 
section of an atom in an atomic state n into a state n' is 
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where the integral is over the natural line width and where the oscillator strength f nn i satisfies 
the sum rule, J2 n ' fnn' — Z, with Z being the number of electrons in the system [13]. If the 
typical energy of the interstellar photons is e ~ 1 eV, and if their energy is boosted in 
the flow's rest frame by a Lorentz factor T ~ 10 3 , then their typical photoabsorption cross 
section by Fe atoms is 

^ fa mn 'hdu = aJ^hcZ ^ iQ _ 18cm ^ 

Te 2m e c Te v 1 

This cross section is larger by about seven orders of magnitude than the typical cross sec- 
tions for production of 7 rays in the interstellar medium by inverse Compton scattering, 
bremsstrhalung and synchrotron emission by electrons, or by n° production by hadrons fol- 
lowed by 7T° — > 27 decay. Thus, if the ionized or excited atoms recombine/decay radiatively 
fast enough, then the typical energy of GRBs (burst size) from NS mergers/AIC in GC is 
given approximately by 

£ 7 ss aN 1 (e x /m A )E K « 3 x 10 51 erg, (6) 

where uia ~ 56m p . Note that the burst size is proportional to the total kinetic energy E K of 
the relativistic flow, and to the column density of the radiation field in the GC along the line 
of sight to the explosion, provided that the flow is not strongly attenuated by the radiation 
field. (If NS merger/AIC occurs near a bright active galactic nucleus, where L ~ 10 12 L Q and 
R ~ 10 17 cm, i.e., iV 7 ~ 10 29 cm~ 2 , then a large enough fraction of the kinetic energy of the 
relativistic flow can be converted into a GRB even by inverse Compton scattering [14,15]. 
However, such a GRB will be extremely short and structureless). 

Let us show that the relativistic atoms are only partially ionized. The recombination 
rate (radiative electron capture) of hydrogen-like atoms into the ground and excited states 
is given approximately by [13] 

r w 4 x 10~ 13 Z 2 T~y /2 cm 3 s-\ (7) 
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where T e y is the temperature in eV in the flow rest frame. The ionization/excitation rate 
must adjust itself to the recombination/decay rate. Consequently, in the rest frame of the 
flow 

— 1/2 

n eU T eV ' Pa cTn^a m 9n 75 /e e y, (8) 

where n e = n e u x 10 11 , n 7 = n 7 5 x 10 5 and r = 1^ x 10 3 . Moreover, in its rest frame the 
relativistic flow expands against the external radiation field until its internal pressure equals 
the external pressure, i.e., 

n e kT ~ Y 2 n 1 e or, n eU T eV ~ Y\n l5 e e v- (9) 

Consequently, T eV m [r 3 e e y/3] 4 / 3 and n e u ^ 9 2 ^ 3 T^ 3 jej^n^. These rather low temper- 
atures and high densities of the relativistic flows justify our initial assumption that the 
inner electronic shells of the relativistic Fe atoms are not ionized and the radiative de- 
cay/recombination is fast enough. 

We have constructed a numerical Monte Carlo code which simulates the production of 
cosmological GRBs by highly relativistic flows in dense stellar regions [16]. The code em- 
ploys the quantum mechanical cross sections for the relevant photo excitation/ionization 
and radiative recombination/decay processes. The core of the MW galaxy [12] was used for 
modeling the stellar environment (density of stars, stellar luminosities and stellar temper- 
atures) in a typical GC. Initial distributions of Lorentz factors which are consistent with 
theoretical considerations and observations have been used. For instance, shock accelera- 
tion usually produces an approximate broken power-law spectrum of Lorentz factors, i.e., 
dn A /dT ~ (T/T m )-P with 1.5 < p < 2.5 for r > T m and p < for T < T m . Using our GRB 
simulation code [16] we have found that (a) the main properties of the simulated GRBs are 
not sensitive to fine details and (b) the calculated light curves and spectral behavior of sim- 
ulated GRBs reproduce remarkably well those observed in GRBs [1, 17-19]. In particular, 
the simulated GRBs look indistinguishable from the observed GRBs. This is demonstrated 
in Fig. 2 which compares a simulated GRB and a GRB from the BATSE IB catalog [18] 
brought in [1] as an example of a typical complex GRB. Note that the temporal power spec- 
tra of the observed and the simulated GRB light curves have the same universal power-law 
behavior, P(w) = \J L(t)exp(iwt)dt\ 2 ~ w~ 2 [15]. This universal power-law behavior can 
be derived analytically [15,16] from our model. Here we summarize briefly approximate an- 
alytical derivations [16] of the other main observed properties of GRBs [1,17-19]. For the 
sake of simplicity we neglect here general relativistic effects (e.g., time dilation and energy 
redshift) and here that the explosions are spherical symmetric and occur at the center of the 
GCs, that the stars within a GC have a uniform spatial distribution, the same luminosity 
and the same effective surface temperature, and that the energy flux in the relativistic flow, 
E 2 driA/dE, is peaked around a Lorentz factor T. Then our model predicts that: 

1. GRB light curves are composed of a smooth background plus strong and weak pulses. 
Strong pulses are produced when the relativistic flow passes near stars within the "beaming 
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cone", i.e., stars at an angle 9* < 1/T relative to the line of sight from the explosion to 
the observer (see eq. 2). When a star is actually a multiple star system (binary, triplet, 
etc) the pulse becomes a multipeak pulse with very short separation in time between the 
peaks (spikes). Weak pulses are produced by boosting star-light from stars near the beaming 
cone. The smooth background is produced by boosting the background light in the beaming 
cone from all the other stars in the GC. Although the photoexcitation of partially ionized 
atoms produces line emission, the line emission is Doppler broadened into a continuum by 
the continuous distribution of the Lorentz factors of the atoms in the flow. Thus, the model 
predicts a continuous energy spectrum that at any moment has an approximate broken 
power-law form (see point 8) which depends on the energy spectrum of the relativistic atoms 
in the flow and on the spectral properties of the star-light which they boost. The ionization 
state of the atoms cuts off emission in the observer frame below E min ~ T m l ~ 5 — 25 keV, 
where I ~ few x 100 eV is the ionization potential of the last bound electrons in the partially 
ionized atoms [16]. It explains why GRBs are not accompanied by detectable X-ray or 
optical-light emission. 

2. The duration of a GRB reflects the spread in arrival times at the observer of gamma 
rays produced within the beaming cone (we neglect the short formation time of the relativistic 
flow because the dynamical time for NS merger/AIC is much shorter, typically 1ms). It is 
given approximately by 

T ~ R/2cT 2 . (10) 

Hence, relativistic flows with F ~ 10 3 in MW-like GC (R ~ 0.5 pc) produce GRBs that last 
typically 25 seconds. Thus, approximately a one light-year path in a GC is contracted into 
a ten seconds 7 ray picture in the observer frame. 

3. The properties of the pulses from single stars depend on their locations, luminosities 
and spectra, and on the distribution of Lorentz factors in the flow. A pulse from a star at a 
distance and an angle begins at a time tj ~ D*6l/2c after the beginning of the GRB 
(t = 0). Its duration (while its intensity is twice the background) is given approximately by 

T p ~ DM0* « b/cT ~ UR/LcF, (11) 

where b is the impact parameter from the star at which the star's light dominates the photon 
column density of the GC. For main sequence stars with 10 _1 < L*/L Q < 10 3 ,, T = 10 3 
and MW-like GC we obtain 5 ms < T p < 50 s. For solar like stars T p ~ 50 ms. Multiple 
stars yield multipeak pulses. For instance, the time difference between the two pulses from 
a binary star is given by 

T b w DM0*/c « d p /Tc, (12) 

where d p is the distance between the binary stars projected on the plane perpendicular to 
the flow. Thus, T b < 500(c?A.u./r3) ms, where c^a.u. is the binary separation in astronomical 
units. Since a large fraction of the stars are in close binaries, triplets, etc, a large fraction of 
the pulses have a multipeak structure. 
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4. In a GRB from a GC with N* stars uniformly distributed, the average number of stars 
within the beaming cone is n p ~ N^/AT 2 . Thus, for a MW-like GC and T ~ 10 3 , we expect, 
on average, n p ~ 3 strong pulses. 

5. The average rate of strong pulses in a GRB for a MW-like GC is dn p /dt ~ cN*/2R ~ 
0.1 s~ l , independent of time. The average time-spacing between the strong pulses, AT, is 
given by, 

AT = T/n p ~ 2R/cN, ~ 10s. (13) 

However, the time-spacing between successive pulses is predicted to fluctuate considerably 
around this average time-spacing. 

6. If the spectrum of T in the flow has a power-law form, driA/dY ~ r~ p , then the energy 
spectrum of a strong pulse is given approximately by [16] 

dnJdE ~ aiT^dE/dT)-^-^ /{T9* + 1/2) (14) 

where a ~ r^ 1 "^. For T > 1/6* E ~ 2re x , 5 < 0.5 one has dn^/dE ~ E^ p+1+S \ For 
r < 1/0*, one has [16] E ~ er 2 , 5 ~ 0, and then dn 7 /d£ ~ #-(p+i+<5)/2_ Hence the S p ectrum 
of a pulse has a broken power-law form. The break occurs when 2T6* ~ 1, i.e., around an 
energy Ef, « 2Fex ~ tx/9* ■ The power-index changes by ~ (p + 1 + 5)/2) ~ 1.5 from 
well below the break to well above the break. Since tj ~ D*9 2 /2c and ~ -R, on average, 
E b ~ y/yJU for small tj and changes to E^ ~ 1/tj for large U. (For a GC with a uniform 
stellar distribution, the probability of is proportional to D 2 and consequently most of 
the stars have & R.) 

7. The relative arrival time at a star of atoms with a Lorentz factor T is given by 
t' = t — ti ~ -D^/2cr -2 . Such atoms boost the star-light to an energy e 7 ~ 2exF. Therefore, 
the peak energy, E p = max E 2 (dn 1 /dE), decreases during a pulse approximately as 

£ p ~ + 5Q a ; 1/2 < a < 1, (15) 

where St* ~ 1 — 100 ms is an added time broadening due to the finite time of the explosion and 
the finite size of the region around the star where the bulk of the photo excitation/ decay takes 
place. E p is maximal right after the pulse begins and decreases monotonically afterwards 
while the photon flux, dn^/dE, peaks at a later time which depends on the distribution of 
Lorentz factors in the flow. Power-law spectra yield photon fluxes during pulses that, on 
average, are time-asymmetric (fast rise and a slower decay) with longer pulses being more 
asymmetric. Pulses are predicted to be narrower and their peak luminosities shifted closer 
to the beginning of the pulse when viewed in higher energy bands (larger Lorentz factors) 
as demonstrated in Fig. 3. 

8. The durations of GRBs have a bimodal distribution which is a trivial consequence of 
their multipulse nature and the fact that, on average, T p <^ AT independent of T: Multipulse 
GRBs have durations which are equal approximately to the sum of the time-spacing between 
their pulses. Consequently, T ~ SAT ^> T p . This produces a bimodal distribution which 
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peaks around T p ~ 0.3 ms for single pulse GRBs and around T ~ 25 s for multipulse GRBs. 
This is demonstrated in fig. 4 for simulated explosions in a MW-like GC. 

9. The photon column densities of MW-like GCs are not large enough to attenuate the 
relativistic flows. When a relativistic flow emerges from a GC it continues to boost star-light 
into gamma rays outside the GC and to emit synchrotron (radio) radiation when it traverses 
interstellar magnetic fields. The extended low level 7 ray emission, which includes much 
broadened and weaker stellar pulses, may last for hours. It is, however, below the current 
detection sensitivity of GRB detectors. The radio emission which lasts for thousands of years 
may be detectable [16]. 

10. The relativistic flows may also collide occasionally with interstellar gas (stellar winds, 
planetary nebulae, etc) or with a molecular cloud with a sizeable column density in/near 
the GC. Typical clouds, have R ~ 5 pc and M ~ 1O 4 M , yielding proton column densities 
of N p ~ 10 22 cm -2 . The total inelastic high energy cross section of iron nuclei on protons 
is a m 10 _24 cm -2 . Thus, nuclear collisions of the flow with a molecular cloud near the 
explosion will produce pions, and consequently, a burst of ~ 10 51 erg multi GeV 7 rays 
through 7T° — > 27 decays (and neutrinos through — > jiu^ and /i — > eu e i>^ decays) with 
a power-law spectrum dn 1 /dE ~ E~ p . For p ~ 2, production of ~ 20 GeV 7 rays comes 
from 7r°'s produced mainly by nuclei with T ~ 200. Thus tt° produced 7 rays of 20 GeV 
are delayed, typically, by D/2cT 2 ~ 2h, as was observed in the case of the 17 February 1994 
GRB [20]. 

In conclusion, we have described a simple mechanism by which NS mergers/accretion 
induced collapse in dense galactic cores produce cosmological GRBs. The remarkable success 
of the model in reproducing all the main observed temporal and spectral properties of GRBs 
[16] strongly suggests that GRBs are 7 ray tomography pictures of dense cores of distant 
galaxies. The proposed mechanism may play an important role also in other astrophysical 
gamma ray sources such as AGN, pulsars and other cosmic accelerators. 
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Figure 1: A schematic drawing illustrating the formation of a GRB by a highly relativistic 
spherical flow in a dense stellar region. Most of the observed 7 rays are produced by the 
radiative decay of photoexcited atoms near stars within a cone with an opening angle 6 ~ 1/T 
along the direction to the observer. 
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Figure 2: (a) A light curve of a simulated GRB and its temporal power spectrum. The 
straight line represents the w~ 2 power law dependence, (b) The light curve of GRB 920110 
from the BATSE IB catalog [18] and its temporal power spectrum. 
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Figure 4: The duration distribution of simulated GRBs from a MW-like GC [16]. 
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